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Abstract 

Atoms trapped in microcavities and interacting through the exchange of virtual photons can 
model an anisotropic Heisenberg spin-1/2 lattice. We do the quantum field theoretical study 
of such a system using the Abelian bosonization method followed by the renormalization group 
analysis. We present interesting physics due to the presence of exchange anisotropy. An infinite 
order Kosterliz-Thouless-Berezinskii transition is replaced by second order XY transition even an 
infinitesimal a small anisotropy in exchange coupling is introduced. We predict a quantum phase 
transition between Mott insulating and photonic superfluid phase due to detuning between the 
cavity and laser frequency, a large detuning favours the photonic superfluid phase. We also do 
the analysis of Jaynes and Cumming Hamiltonian to support results of quantum field theoretical 
study. 
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I. INTRODUCTION 



The physics of strongly correlated system is interesting in its own right and manifests 



in different branches o 
natural oxide materials 



physics. Some of the important correlated physics appears in 
1] and some of them are in engineered materials, like the correlated 
ohysics in Josephson junction array |2j, Bose-Einstein condensation and optical lattice 
3), 4|. Therefore one can raise the question what is the further source of correlated 
physics in the state of enginnering ?. The recent experimental success in engineering 
strong interaction between the photons and atoms in high quality micro- cavities opens up 



;he possibility to use 
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ight matter system as quantum simulators for many body physics 
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Here we would like to discuss the basic physics of micro-optical cavity very briefly. A micro- 
cavity can be created in a photonic band gap material by producing a localised defect in the 
structure of the crystal, in such a way that light of a particular frequency can not propage 
outside the defect area. Large arrays of such micro-cavities have been produce. Photon hop- 
ping between neighboring cavities has been observed in the microwave and optical domains. 
Many body Hamiltonians can be created and probed in coupled cavity arrays. There atoms 
are used for detection and also for the generation of interaction between photons in the same 
cavity. As the distance between the adjacent cavities is considerably larger than the optical 
wave length of the resonant mode, individual cavities can be addressed. These artificial 
system could act as a quantum simulator. In this optical cavities system we study the 
different quantum phase of polariton ( a combined excitations of atom-photon interactions. 
) by using spin models that conserve the total number of excitations. Best of our knowledge, 
at first we have done the explcit quantum field theoretical calculations of this type of system. 



At first we would like to discuss the generation of the spin model for such type of systems. 
It has already been discussed in the literature but we mention this process very briefly for 
the sake of completeness. In this description we will follow the Ref. [f| and Ref. 3]. 

Micro-cavities of the photonic crystal are coupled through the exchange of photons. Each 
cavity consists of one atom with three levels in the energy spectrum, two of them are long 
lived and represents two spin states of the system and the other is excited. Externally 
applied laser and cavity modes couple to each atom of the cavity. It may induce the Raman 



2 



transition between these two long lived levels. Under a suitable detuning between the laser 
and the cavity modes virtual photons are created in the cavity which mediate an interaction 
with another atom in a neighboring cavities. One can eliminated the excited states of atomic 
level photon states by choosing the appropriate detuning between the applid laser and cavity 
modes. Then one can achieve only two states per atom in the long lived state and the system 



can be described by a spin-1/2 Hamiltonian 



Fig. 1 shows the schematic phase diagram of our description to generate the xy spin 
interaction of the system. The Hamiltonian of the system consists of three parts: 

H = H A + H B + H c (1) 

Hamiltonians are following 

N 

H A = ^uJe\ej >< ej\ + u ab \bj >< bj\ (2) 

3=1 

j is the cavity index. u a b and u e is the energy of the state \b > and the excited state 
respectively. The energy level of state \a > is set zero. The following Hamiltonian describes 
photons in the cavity, 

N N 

H c = uj c J2 a j ta j + ^^(a/ij+i + M. ( 3 ) 

3=1 3=1 

where af creates a photon in cavity j, uoq is the energy of photons and Jc is the tunneling 
rate of photons between neighboring cavities. Interaction between the atoms and photons 
and also by the driving lasers are described 

N n 

Hac = EK -fe-^ + ^ a a i )|e i >< aj\ + h.c] + [a^b}. (4) 

3=1 1 

Here g a and g>& are the couplings of the cavity mode for the transition from energy states \a > 
and 1 6 > to the excited state. Q a and Qf, are the Rabi frequencies of lasers with frequencies 
u a and u h respectively. 

They have derived an effective spin model by considering the following physical processes: 
A virtual process regarding emission and absorption of photons between the two stable states 
of neghiboring cavities. The resulting effective Hamiltonian is 

N N 

H xy = Yj Ba 3 Z + J2( J l a ^ a 3+l~ + J 2VfCTj+l~ + h.c) (5) 
3=1 3=1 



When J 2 is real then this Hamiltonian reduces to the XY model. Where of = \bj >< 

bj\ ~ \ a j >< a j\i a j + = \bj >< a j\, °f = \ a j >< bj\ 



N 



Hxy = ^2B(ai Z + J x ai x a i+ i x + J y ai y a i+ i y ). (6) 
With J x = (J 1 + J 2 )/2 and J y = (J x - J 2 )/2. 

Here we discuss very briefly about an effective ofo* in such a system. Authors of 
Ref. have proposed the same atomic level configuration but having only one laser of 
frequency uj that mediates atom-atom coupling through virtual photons. Another laser 
field with frequency v is used to tune the effective magnetic field. They described the one- 
dimensional case. In this case the Hamiltonian Hac will change but the Hamiltonians Ha 
and He will not. 

N £1 A 

Hac = EK^T e " W + T^-^a^e, >< a,| + h.c] + [a <- 6]. (7) 

Here, fl a and fib are the Rabi frequencies of the driving laser with frequency u on transition 
I a >— > |e > 1 6 >— > |e >, whereas A a and Aj, driving laser with frequency z/ on transition 
I a >— > \e > \b >^ \e >. They have eliminated adiabatically the excited atomic levels and 
photons by considering the interaction picture with respect to Ho = Ha + Ho- They 
have considered the detuning parameter in such a way that the Raman transitions between 
two level supressed. They have also chosen the parameter in such a way that the dominant 
two-photon processes are those that involve one laser photon and one cavity photon each 
but the atom does no transition between levels a and b. Whenever two atoms exchange a 
virtual photon then both of them experience a Stark shift plays the role of an efective o z o z 
interaction. Then the effective Hamiltonian reduce to 

N 

H zz = + J s a j 'a j+1 z ) (8) 

Analytical expressions for B, J 1; J 2 and J z has given in Ref. These two parameters can 
be tuned independently by varying the laser frequencies. They have obtained an effective 
model by combining Hamiltonians H xy and H zz by using Suzuki- Trotter formalism. The 
effective Hamiltonian simulated by this procedure is 

TV 

(9) 

3=1 a=x,y,z 



where B to t = B + B. It has been shown in Ref. [7| that J y is less than than J x . It is clear 
from analytical expressions for J x and J y that the magnitudes of J\ and J 2 are different. 
This result of numerical simulations trigger us to define a model, which has given below to 
study the quantum phases of this system and also the transition among them and at the 
same time this subject is in the state of art of engineering. 



II. RENORMALIZATION GROUP STUDY OF MODEL HAMILTONIAN 

We consider the anisotropic Heisenberg spin-1/2 Hamiltonian on a one dimensional lat- 
tice. The XYZ Heisenberg Hamiltonians is defined as: 

Hxyz = £ i + + (l-a)SlSl +l 

n 

+ A S z n S z n+1 + hS z n ], (10) 

where are the spin-1/2 operators. We assume that the XY anisotropy a and the zz 
coupling A satisfy — 1 < A < 1, and < a < 1 and magnetic field strength h > 0. The 
Hamiltonian H X yz is invariant under the transformation — ► —S%, — > —S^, S* — > —S^, 
actually it is a Z 2 symmetry. For finite h, Z 2 symmetry is absent when — > —S^. Here 
h ~ B tot , A = J z , Ji = 1 and J 2 = a. 

Spin operators can be recasted in terms of spinless fermions through Jordan- Wigner trans- 
ormation and then finally one can express the spinless fermions in terms of bosonic fields 



19]. We recast the spinless fermions operators in terms of field operators by this relation 
19] . 

1>(x) = [e lkpx Mx) + e~ lkFX Mx)] ( n ) 

where i[>r(x) and describe the second-quantized fields of right- and the left- moving 

fermions respectively, hp is Fermi wave vector. Therefore, one can study the effect of gate 
voltage through arbitrary hp. We would like to express the fermionic fields in terms of 
bosonic field by the relation 

^ r ( x ) = g -< - ^ ( 12 ) 

y2ixa 

r is denoting the chirality of the fermionic fields, right (1) or left movers (-1). The operators 
U r commutes with the bosonic field. U r of different species commute and U r of the same 
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species anti-commute. field corresponds to the quantum fluctuations (bosonic) of spin and 
6 is the dual field of <p. They are related by the relations 4>r = 6 — (p and 4>l = 6 + (p. 
Hamiltonian H is non-interacting part of Hxyz- 



Ho = - j dx [ (d x ey + (d x <i>y ] , (13) 

where v is the velocity of the low-energy excitations. It is one of the Luttinger liquid 
parameter, and the other is K, which is related to A by Q 

K = a + 2 3 in--(A) ' (14) 

K takes the values 1 and 1/2 for A = (free field), and A = 1 (isotropic anti-ferromagnet), 
respectively. The relation between K and A is not preserved under the renormalization, so 
this relation is only correct for initial Hamiltonian. The analytical form of the spin operators 
in terms of the bosonic fields are 

~7T 



SI = [c 2 cos(2V^0 + (-l) n c 3 ] 



SI = -[c 2 cos(2vW0) + (-1)%] sin(y-0), 

S z n = \f^d x( f> + (-l) n c lC os(2V^K<f>) , (15) 



where qs are constants as given in Ref. [20(. The Hamiltonian Hxyz i n terms of bosonic 
fields is the following, 



Hxyz = H + a J cos(2^J^9(x))dx + A J cos(4V 7iK(p(x))dx 

—h / d x (p(x)dx (16) 



One can get the Hxy Hamiltonian by simply putting A = in the above Hamiltonian. 
In this derivation, different powers of coefficients q have been absorbed in the definition 
of a, h and A. The integration of the oscillatory terms in the Hamiltonian yield negligible 
small ontributions, the origin of the oscillatory terms are the spin operators. So it's a 
resonably good approximation to keep only the non-oscillatory terms in the Hamiltonian 
. The Gaussian scaling dimension of these coupling terms, a, A are 1/K, AK respectively. 
The third term (A) of the Hamiltonian tends to order the system into density wave phase, 
whereas the second term (a) of the Hamiltonian favours the staggered order in the xy 
plane. Two sine-Gordon couplings terms are from two dual fields. Therefore the model 



Hamiltonian consists of two competing interactions between the ordered phase and the XY 
order. This Hamiltonian contains two strongly relevant and mutually nonlocal perturbation 
over the Gaussian (critical) theory. In such situation strong coupling fixed point is usually 
determined by the most relevant perturbation whose amplitude grows up according to its 
Gaussian scaling dimensions and it is not much affected by the less relevant coupling terms. 
However this is not the general rule, if the two operators exclude each other, i.e., if the field 
configurations which minimize one perturbation term do not minimize the other. In this 
case interplay between the two competing relevant operators can produce a novel quantum 
phase transition through a critical point or a critical line. Therefore we would like to study 
the RG euation to interpret the quantum phases of the system. 

We will now study how the parameters a, A and K flow under RG. The operators in Eq. 
17 are related to each other through the operator product expansion; the RG equations for 
their coefficients will therefore be coupled to each other. Here we derive the RG equations 
by using perturbative renormalization group approach scheme. We use operator product 
expansion to derive these RG equations which is independent of boundary condition 24]. 
RG equations themselves have been established in a perturbative expansion in coupling 



constant (g(l)), they cease to be valid beyond the certain length scale, where g(l) ~ 1 19]. 
The RG equations for the coefficients of Hamiltonian Hxyz are 

These RG equations have trivial (a* = = A*) fixed points for any arbitrary K. Apart 
from that these RG equations have also two non-trivial fixed lines, a = A and a = —A for 
K = 1/2. In our study, there are critical surfaces on which the system flows onto the non- 
trivial fixed lines (a = ± A) . A density wave states can be characterized when K — > or the 
staggered ordered when (K — > oo). Note that the transition occuring on them are second 
order. Infinitesimal amount of anisotropy will change the situation drastically, a gapless 
phase in absence of anisotropy will change to the gapped phase in presence of anisotropy. 
Since this gapped excitation is not directly related to magnetization, therefore it will not 
favour to create the plateau phase. When the system is in the plateau phase the transition 
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driven by the magnetic field is always of z = 2 ( z is the dynamical critical exponent) and 
thus the plateau shows a square-root behaviour of magnetization. When a is increasing, a 
second order transition drives system to the in plane XY ordered phase , whose exponents 
depend on initial couplings and hence are nonuniversal. In absence of planar anisotropy 

£osterliz- Thouless-Brezinskii (KTB). When the inplane 



the transition to plateau state is 
anisotropy is present, then z = 1 [2JJ. Please see Refs. |22|, |23[ for detail understanding of 
this subject. A magnetic field larger than the relevant gap of the system drives the system 
to a gapless phase. This transition is from commensurate phase to incommensurate phase 
transition. 

We have seen the analytical expression of B tot from Ref. [7] that the total magnetic field 
increase for the larger values of detuning, therefore larger detuning drives the system from 
gapped (Mott-insulating) state to gapless superfluid state. 

We now discuss, how the effective repulsion will decrease as we increase the detuning 
between the atomic and laser frequency. It can be explained starting from the Jaynes- 



Cummings Hamiltonian 25l. l26| . Janes- Cummings Hamiltonian for a single atom is 



Hjc = uca^a + uj \e >< e\ + X(a^\g >< e\ + a\e >< g\) (18) 

toe and ujq are the frequency of the resonant mode of the cavity and of the atomic transition, 
respectively. A is the Jaynes-Cumming coupling between the cavity mode and the two level 
system, (a) is the creation (annihilation) operator of a photon inside a cavity. \g > and 
|e > is respectively the ground state and excited states of the two level system respectively. 
When we consider large values of photon and atomic frequencies compare to atom-photon 
coupling A, the number of excitations is conserved for this Hamiltonian. Suppose we consider 
fixed numbers of excitation, n. The energy eigenvalues for n excitations is 



A A 2 
E n =nu c + -± y nA 2 + — (19) 

Here A = (uj — uj c ) and n > 1. Now we consider an array of cavities, the basic Hamiltonian 
for each cavity is the same as that of Eq.19 . Here we consider the system with one excitation 
of energy E\ in each cavity and the lowest energy of two excitation in each cavity is E 2 . 
Therefore to create one additional excitation in each cavity requires energy 



A 2 A 2 A 

E 2 -2 El =2^ + — - v/2A 2 + 
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Which one may consider as an effective on-site repulsion because it measures the difference 
between the energy of two and single excitation (polariton in each cavity). This effective 
repulsion decreases as we increase the detune factor. Therefore we can conclude that for 
A = 0, double occupation never occurs, indicating a Mott insulating behaviour. When A is 
much larger than the coupling A then the occupation number larger than one occurs as to be 
expected for a photonic superffuid regime. In our quantum field theoretical calculations, we 
have also predicted that the large detuning drives the system from gapped Mott insulating 
phase to the gapless superffuid phase. 

CONCLUSIONS 

At first we have done the quantum field theoretical analysis of an effective spin model 
in coupled optical cavity arrays. We have predicted two quantum phases, Mott insulator 
and photonic superfiudity. Anisotropy in the exchange interaction has also created a gapped 
phase. An infinite order KTB transition has been replaced by the second order XY transition. 
The rigorous quantum field theoretical derivation of this manuscript is absent in all previous 
studies and also we provide physical explanation of the transition process based on Jaynes- 
Cummings Hamiltonian. 
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